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/  ABSTRACT 

Due  to  the  lack  of  longitudinal  stiffening  along 


submarine  hulls,  they  are  inherently  vulnerable  to 
underwater  explosions  directly  off  the  bow  or  stern. 
Accordion-like  deformations  of  the  hull  are  set  into  motion 
which  could  cause  dynamic  amplification  in  the  transient 
response  of  attached  substructures.  In  underwater  shock 
acceptance  tests  of  internal  equipment,  this  interaction  is 
created  by  exploding  a  charge  in  a  fore  and  aft  configura¬ 
tion  with  the  submerged  shock  test  vehicle  (SSTV) .  With  the 
increasing  availability  of  large  computers  and  the  rapid 
development  of  numerical  methods,  several  computer  codes 
have  been  written  to  predict  equipment  response  to 

underwater  shocks.  Using  the  ELSHOK  (ELASTIC  SHOCK)  code, 

(  I  j  rs 

this  investigation  studied  the  effect  of  hull/substructure 
interaction  on  stiffened  shell  response  at  resonance 
following  an  end-on  load.  The  transient  response  of  the 
coupled  shell/substructure  system  from  tapered  and 
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I.  INTRODUCTION 

A.  BACKGROUND 

Submarines  are  inherently  vulnerable  to  explosions 
detonated  directly  off  the  bow  or  stern.  While  they  have 
adequate  protection  against  side  loads  with  transverse 
bulkheads  and  frames,  there  is  little  longitudinal 
protection  along  the  hull.  As  a  result  of  this,  they  are 
susceptible  to  accordion-like  deformations  of  the  hull  upon 
impact  from  end-on  loads.  It  is  possible  that  the  frequency 
of  the  hull  motions  could  excite  internal  equipment  into 
resonance  thereby  causing  unacceptable  damage.  This 
investigation  examines  the  coupled  elastic  response  of  a 
shell/substructure  system  when  subjected  to  an  end-on 
underwater  explosion  (UNDEX)  as  the  mass  and  stiffness  of 
the  internal  substructure  is  varied. 

The  current  specifications  for  conducting  underwater 
shock  tests  on  submar ine- installed  equipment  are  contained 
in  MIL-S-901D  [Ref.  1],  This  document  specifies  the 
explosive  charge  weight  and  geometry  of  the  test,  as  well  as 
the  mounting  and  orientation  of  the  equipment  being  tested 
within  the  submerged  shock  test  vehicle  (SSTV) .  The  SSTV  is 
nothing  more  than  a  ring-stiffened  cylindrical  shell  with 
circular  endplates  designed  to  simulate  the  hull  motions 


created  by  an  underwater  shock.  Recognizing  the  potential 
for  unacceptable  damage  to  hull  mounted  equipment  from  an 
end-on  load,  the  specification  requires  the  first  of  four 
UNDEX's  to  be  performed  in  a  fore  and  aft  configuration. 

The  suitability  of  an  equipment  design  or  installation  is 
evaluated  according  to  its  ability  to  function  as  intended 
during  and  after  each  shock  impulse.  Equipment  tests  of 
this  nature  are  quite  expensive  and  require  a  great  deal  of 
preparation.  It  is  desirable  for  the  designer  to  have  some 
idea  of  what  the  transient  response  of  the  hull  and 
equipment  will  be  prior  to  the  actual  UNDEX. 

The  development  of  analytical  methods  and  computer  codes 
for  the  analysis  of  the  reaction  of  submerged  structures  to 
underwater  explosions  has  enabled  the  engineer  to  simulate 
the  shock  response  of  submarines  and  equipment  with 
increasing  accuracy.  Finite-element/finite  difference 
methods  allow  for  the  analysis  of  structure  and  fluid 
responses,  and  the  work  of  Geers  [Ref.  2]  provides  a  means 
to  incorporate  fluid-structure  interaction  effects.  Several 
computer  codes  currently  exist  which  utilize  these 
principles  to  analyze  installed  equipment  response  to  shock 
waves.  One  of  them  is  the  ELSHOK  (ELASTIC  SHOCK)  code 
developed  by  Weidlinger  Associates,  under  the  sponsorship  of 
the  Defense  Nuclear  Agency  (DNA)  and  the  Office  of  Naval 
Research  (ONR) .  This  code  was  developed  to  investigate 
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modern  submarine  underwater  explosive  shock  response  in 
conjunction  with  a  testing  program  using  small  to  large 
scale  models  and  shaped  or  tapered  explosive  charges.  The 
accuracy  of  ELSHOK  has  been  validated  in  several  highly 
controlled  tests,  notably  the  1983  low  level  explosive  test 
of  an  SSN  668  class  submarine.  The  code  was  made  available 
to  the  Naval  Postgraduate  School  by  the  DMA  with  support 
from  Weidlinger  Associates.  It  was  first  used  at  the  school 
by  LT  Mark  Welch,  USN,  in  his  study  comparing  the  shock 
response  predictions  of  ELSHOK  with  those  obtained  using  the 
Dynamic  Design  Analysis  Method  (DDAM)  [Ref.  3]. 

B.  PURPOSE  FOR  THIS  INVESTIGATION 

An  explosive  detonated  underwater  exerts  a  great  amount 
of  pressure  on  the  surrounding  fluid.  The  water  is  actually 
compressed  by  the  force  of  the  explosion  causing  a  shock 
wave  to  form  at  the  point  of  detonation  which  propagates  in 
a  roughly  spherical  shape  at  the  speed  of  sound  in  water. 
When  the  pressure  pulse  impinges  on  a  submarine,  dynamic 
responses  result  in  the  hull  from  the  fluid-structure 
interaction.  The  submarine  experiences  translational  motion 
away  from  the  point  of  impact.  The  velocity  of  this  motion 
is  of  particular  interest  to  a  designer  of  submarines  or 
weapons . 

In  an  explosion  with  a  side-on  aspect  with  the 


submarine,  the  shock  wave  hits  the  hull  in  the  middle  before 


It  hits  the  ends  due  to  its  spherical  shape.  This  causes 
the  middle  to  bow  away  from  the  point  of  impact  while  the 
ends  remain  initially  fixed.  The  resulting  motion  is  called 
"whipping"  as  the  hull  bends  back  and  forth  about  the 
middle.  With  an  explosion  off  the  bow  or  stern,  the  shock 
wave  hits  the  submarine  as  a  plane  wave  due  to  the  small 
cross-sectional  area  at  the  point  of  impact.  Since  the 
force  of  the  shock  wave  is  axisymmetric  about  the  centerline 
of  the  submarine,  no  whipping  motion  is  induced  in  the  hull. 
The  inertial  forces  are  transmitted  down  the  hull  until  they 
reach  the  stern.  The  stern  reacts  violently  to  the  combined 
inertial  forces  resulting  from  the  explosion  as  well  as  the 
reflected  pressure  pulses  from  the  water  displaced  by  the 
translational  motion  of  the  submarine.  After  the  initial 
transient  response,  the  hull  settles  into  an  "accordion 
mode"  about  the  middle  where  the  motion  of  the  bow  opposes 
that  of  the  stern. 

The  purpose  of  this  investigation  is  to  examine  the 
shell/substructure  interaction  when  subjected  to  UNDEX  end- 
on  loading  using  the  ELSHOK  code.  The  hull  used  in  this 
study  is  a  ring-stiffened  cylindrical  shell  with  endplates 
similar  in  shape  to  those  used  as  SSTV's.  The  internal 
equipment  or  substructure  is  a  diaphragm  attached  to  the 
shell  whose  thickness  is  increased  in  order  to  study  how  the 
mass  ratio  between  substructure  and  shell  affects  the 


coupled  transient  response.  The  dynamic  amplification  of 
the  response  of  the  substructure  at  resonance  is  examined 
for  both  taper  and  conventional  charges  of  equivalent 
impulse. 


EXPLOSIVE  SHOCK RESPONSE  USING  ELSHOK 


A.  GENERAL  PRINCIPLES  OF  OPERATION 

The  ELSHOK  computer  code  consists  of  a  family  of 
programs  developed  to  calculate  the  transient  response  of  a 
submerged,  ring-stiffened  shell  of  revolution  of  finite 
length  to  an  underwater  shock  wave  emanating  from  an 
explosive  source  located  at  an  arbitrary  point  away  from  the 
structure  [Ref.  4].  The  shell  is  assumed  to  be  linearly 
elastic,  with  or  without  internal  substructures,  and  the 
surrounding  fluid  is  treated  as  an  infinite  acoustic  medium. 
Component  modal  analysis  is  employed  in  all  phases  of  the 
calculations.  The  complete  structural  system  is  considered 
to  consist  of  the  ring-stiffened  shell  and  any  attached 
substructures.  The  vibration  modes  of  each  component  are 
calculated  separately,  and  the  equations  of  motion  for  the 
entire  system  are  obtained  by  enforcing  compatibility  of 
deformation  at  the  points  of  attachment.  The  free-free 
modes  of  the  empty  ring-stiffened  shell  and  the  fixed-base 
modes  of  each  individual  substructure  are  coupled  through 
the  use  of  dynamic  boundary  conditions  [Refs.  5,  6].  This 
eliminates  the  need  to  calculate  the  combined  modes  and 
natural  frequencies  of  the  entire  system  as  well  as  the 
requirement  for  a  combined  system  stiffness  matrix. 


The  structure-fluid  interaction  is  approximated  in 
BLSHOK  using  the  Doubly  Asymptotic  Approximation  (DAA) 
methods  of  Geers  [Ref.  2],  expressed  in  terms  of  functions 
which  are  orthogonal  over  the  wet  surface  of  the  submerged 
shell.  By  matching  exact  pressure-velocity  relations  at 
zero  and  infinite  frequencies,  the  elements  of  the  matrices 
in  the  resulting  DAA  are  obtained.  In  transient  problems, 
the  DAA  yields  exact  solutions  at  early  and  late  times 
providing  a  smooth  transition  between  these  two  limits.  By 
accounting  for  the  effects  of  the  fluid  with  quantities 
defined  solely  on  the  wet  surface  of  the  shell,  the  DAA 
essentially  uncouples  the  fluid  field  from  the  structural 
field . 

B.  ORGANIZATION  AND  OPERATION  OF  ELSHOK 

As  stated  previously,  the  ELSHOK  code  contains  a  group 
of  computer  programs  which  are  utilized  in  order  to  obtain 
velocity-time  histories  for  various  locations  in  the  shell- 
substructure  system.  The  programs,  in  order  of  execution, 
are : 

1)  B0S0R4 — structural  analyzer  for  shell  [Ref.  7] 

2)  ACESNID — virtual  mass  processor 

3)  PIFLASH — shell-fluid  processor 

4)  SAPIV — structural  analyzer  for  substructure  [Ref.  8] 

5)  PICRUST — substructure  processor 

6)  USLOB — time  integration  processor 


7)  PUSLOB — plotting  processor 


Figure  1  illustrates  the  general  relationship  between  the 
programs  which  make  up  the  ELSHOK  computer  code.  As  can  be 
seen,  the  operation  is  divided  into  four  phases. 

1 .  Phase  I — Shell  and  Fluid  Analysis 

The  first  step  of  the  ELSHOK  calculation  is  to 
determine  the  in-vacuo  free-free  modes  and  natural  frequen¬ 
cies  of  the  shell.  This  is  accomplished  by  modeling  the 
shell  using  the  B0S0R4  finite  difference  code.  This  code 
applies  to  segmented,  ring-stiffened,  branched  shells  of 
revolution  having  various  meridional  geometries,  wall 
constructions,  and  ring  reinforcements  making  it  well  suited 
for  the  modeling  of  submarines  or  SSTV's.  For  compatibility 
with  ELSHOK,  a  separate  B0S0R4  calculation  must  be  performed 
for  each  circumferential  harmonic  (N)  included  in  the 
analysis.  Examples  of  different  harmonics  in  the  circum¬ 
ferential  distributions  include  N  =  0  (breathing/torsional) 
and  N  =  1  (pure  translation/whipping)  as  shown  in  Figure  2. 

The  second  part  of  Phase  I  is  the  computation  of -the 
virtual  mass  array  using  ACESNID.  This  provides  the 
late-time  contribution  of  the  DAA.  The  virtual  mass  array 
is  determined  from  the  solution,  based  on  simple  sources  of 
a  low-frequency  steady-state  problem  in  which  normal 
displacements  corresponding  to  surface  expansion  functions 
are  applied  to  the  surface  of  revolution  in  the  infinite 


cibcumferential  nooal  pattern 


Figure  2.  Nodal  Patterns  for  Cylindrical  Shells 


fluid  having  the  same  shape  and  size  as  the  net  surface  of 
the  shell.  ACESNID  requires  only  one  execution  for  all 
values  of  N  considered. 

The  execution  of  PIFLASH  completes  Phase  I  by 
combining  the  outputs  from  BOSOR4  and  ACESNID  to  create  a 
"shell-fluid  file"  which  contains  all  the  information 
required  to  describe  the  shell  during  subsequent  operations. 
One  of  the  useful  features  of  ELSHOK  is  the  fact  that  any 
number  of  substructures  can  be  analyzed  for  a  given  shell 
without  having  to  repeat  Phase  I  calculations. 

2 .  Phase  II — Substructure  Analysis 

The  analysis  of  the  substructure  is  performed  using 
the  SAPIV  finite  element  code.  The  SAPIV  code  contains  a 
variety  of  finite  elements  (plates,  beams,  pipes,  etc.) 
which  can  be  used  in  the  modeling  process.  It  is  a  general 
purpose  code  based  upon  a  lumped  mass  formulation  applying 
to  linearly  elastic  structures.  The  PICRUST  code  takes  data 
from  SAPIV  and  reorganizes  it  to  facilitate  the  solution  of 
the  equations  governing  the  transient  response  problem  under 
study.  It  is  during  the  execution  of  PICRUST  that  the 
connectivity  between  the  substructure  and  the  shell  is 
accounted  for  and  influence  coefficients  for  evaluating  the 
forces  developed  at  the  points  of  attachment  are  calculated. 


3.  Phase  III — Submerged  Shock  Response 


shell-fluid  file  from  Phase  I  is  merged  with  the 
substructure  files  from  Phase  II.  OSLOB  is  the  portion  of 
the  ELSHOK  code  which  contains  the  underwater  shock  from  an 
arbitrary  point.  It  allows  pressure-time  history  inputs  for 
taper  charge  modeling  or  an  exponentially  decaying  pressure 
impulse  such  as  those  encountered  with  conventional  charges. 
It  is  of  the  form  [Ref.  9]: 

P(t)  =  Kj^(w^/^/R)  ^  exp(-t/eQ)  (1) 

where : 

P(t)  =  incident  pressure  on  the  shell  (psi) 

=  multiplicative  constant  for  incident  pressure 

K2  =•  spatial  decay  constant  for  incident  pressure 

t  =  time  after  arrival^of  shock  wave  at  point  of 

interest  (msec) 

w  =  weight  of  spherical  charge  (lb.) 

tr 

1/2  1/2  4 

®o  “  (w  /R)  =  time  constant  of  exponential 

decay  (msec) 

®  multiplicative  constant  for  time  constant 
=  spatial  decay  constant  for  time  constant 

R  =  distance  from  the  explosive  to  the  point  of 
interest  (ft) 

Kj^,  ^2!  K^/  and  are  constants  which  depend  on 
explosive  type.  USLOB  employs  a  modified  version  of  the 
Runge-Kutta  integration  method  to  produce  velocity-time 
histories  in  tabular  form  at  user  specified  points  on  the 
shell  and  the  substructure. 


4 .  Phase  IV — Plots  of  Velocity-Time  Histories 
After  the  velocity-time  histories  have  been 
obtained,  PUSLOB  is  utilized  to  display  the  velocity 
responses  in  a  plotted  format  on  a  TEKTRONIX  graphics 
terminal.  These  plots  are  adequate  for  viewing  trends  and 
general  responses,  however  for  formal  display  purposes,  the 
IBM  Versatec  plotter  yields  a  superior  product.  Appendix  A 
contains  a  code  which  was  written  to  convert  the  velocity 
punch-card  files  created  by  PUSLOB  into  a  data  format  which 
is  acceptable  for  Versatic  plotting  using  the  EASYPLOT 


program 


Ill 


MODEL  USED  IN  THE  ANALYSIS 


A.  SHELL  MODEL 

The  shell  used  in  this  analysis  is  a  small-scale  model 
of  a  typical  SSTV  used  for  conducting  underwater  shock  tests 
on  high-impact  shipboard  machinery.  As  can  be  seen  in 
Figure  3,  it  is  a  high  strength  steel,  free-free  ring 
stiffened  cylinder  with  aluminum  endplates.  Appendix  B 
includes  a  sample  B0S0R4  input  for  one  of  the 
circumferential  harmonic  distributions.  Data  can  be  input 
to  B0S0R4  using  a  formatted  input  code  or  by  responding 
interactively  to  a  set  of  prompts.  Each  segment  is  modeled 
separately  which  allows  for  different  materials  and  physical 
properties  to  be  included  in  the  structure.  The  small 
stiffeners  along  the  shell  are  represented  by  an  orthotropic 
approximation  which  in  effect  increases  the  density  of  the 
shell.  The  six  discrete  rings  are  modeled  individually, 
allowing  for  different  size  or  material.  The  torsional 
rigidity  (GJ)  for  each  ring  is  found  using  the  following 
relation; 

GJ  =  EJ/2 (1  +  V)  (2) 

where ; 

E  =  modulus  of  elasticity  (psi) 
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J  =  A  /(40  Ip)  =  torsional  rigidity  constant  [Ref.  10] 

V  =  Poisson's  ratio 

Since  the  shell  is  a  surface  of  revolution,  the 
properties  of  each  ring  cross-section  need  only  be  specified 
at  one  point.  The  end  plates  (segments  one  and  three)  are 
divided  into  ten  equally  spaced  nodes,  and  the  cylindrical 
portion  of  the  shell  (segment  two)  is  divided  into  forty- 
five  nodes  along  the  longitudinal  axis.  As  with  any  finite 
difference  code,  greater  accuracy  will  be  obtained  in  the 
solution  with  higher  numbers  of  nodes,  however,  computa¬ 
tional  time  will  increase.  Two  additional  nodes  are 
automatically  inserted  by  BOSOR4  into  each  segment  in  order 
to  reduce  the  truncation  errors  associated  with  segment 
interfaces  and  to  prevent  spurious  vibration  or  buc)cling 
modes.  Through  finite  difference  techniques,  B0S0R4 
determines  the  in-vacuo  free-free  modes  and  natural 
frequencies  of  the  shell.  Since  this  study  deals  strictly 
with  end-on  loading,  the  problem  is  axisymmetric  about  the 
longitudinal  axis.  Although  only  the  N  =  0  (breathing/ 
torsional)  modes  are  activated,  a  few  N  =  1  (translational) 
models  must  be  calculated  and  retained  due  to  the  fact  that 
the  ELSHOK  computer  code  was  developed  for  side-on  loading 
and,  therefore,  requires  a  translational  input.  Since  there 
is  no  translation,  the  solution  is  not  affected  by  these 
superfluous  modes.  BOSOR4  generates  a  shell  file  which 


gives  the  mode  shapes  for  each  natural  frequency.  All 
breathing  modes  have  been  retained  for  further  study, 
however,  the  torsional  modes  have  been  dropped  from  the 
analysis  since  they  are  not  activated  by  an  end-on  load. 
Table  I  contains  a  listing  of  modes  and  natural  frequencies 
retained  for  input  into  the  PIFLASH  program. 

Figure  4  gives  a  representation  of  a  typical  shell- 
substructure  configuration.  The  two  coordinate  systems  used 
in  an  ELSHOK  analysis  are  illustrated.  X,  Y,  and  Z  refer  to 
the  global  coordinate  system  of  the  shell,  while  x,  y,  and  z 
refer  to  the  local  coordinate  system  of  the  substructure. 

In  both  systems,  the  x-axes  run  longitudinally  down  the 
hull.  The  substructure  z-axis  is  related  to  the  global 
Z-axis  by  an  angle  where  a  refers  to  the  substructure. 

r 

Points  on  the  wet  surface  of  the  shell  are  located  by 
specifying  the  meridional  arc  length  s,  measured  along  the 
reference  surface,  and  the  circumferential  coordinate  angle 
0.  The  lower  case  u,  v,  and  w  are  used  for  local  shell 
displacements.  To  eliminate  torsional  modes  from 
consideration,  all  N  =  0  modes  in  the  shell  file  with 
significant  v-displacements  are  discarded. 

Once  the  shell  files  have  been  obtained  from  B0S0R4,  the 
input  code  for  ACESNID  can  be  created.  This  program 
determines  the  virtual  mass  array  which  produces  the  late¬ 
time  contribution  of  the  DAA,  It  considers  a  cavity  in  an 


TABLE  I 


MODES  AND  NATURAL  FREQUENCIES  OF  SHELL 


infinite  acoustic  fluid  corresponding  to  the  wet  surface  of 
the  shell.  A  user  specified  number  of  surface  expansion 
functions  having  the  property  of  orthogonality  over  the  wet 
surface  of  the  shell  are  obtained  during  the  calculation  of 
the  virtual  or  entrained  mass.  The  cavity  is  divided  into 
bards  which  describe  the  behavior  of  the  fluid.  A 
sufficient  number  of  bands  must  be  provided  to  enable  the 
fluid  response  to  match  the  normal  motion  specified  by  the 
surface  expansion  functions.  In  a  sense,  they  can  be  looked 
upon  as  nodes  in  a  finite  element  model.  In  this  analysis, 
the  shell  is  assumed  to  be  totally  immersed  in  salt  water 
with  a  mass  density  of  9.59684E-05  Ibf-s^/in^.  The 
endplates  are  divided  into  fifteen  bands  each,  and  the 
cylinder  has  fifty-one  bands.  Six  surface  expansion 
functions  have  been  generated  for  the  endplates  with  ten 
along  the  cylinder. 

The  shell  mode  files  from  BOSOR4  and  the  virtual  mass 
file  from  ACESNID  are  combined  in  the  PIFLASH  program  to 
produce  a  shell-fluid  file  for  further  processing.  The 
model  of  the  shell  is  now  complete  and  can  be  used  for  any 
variety  of  substructures  without  any  recalculation. 

B.  SUBSTRUCTUPE  MODEL 

The  idealized  internal  equipment  selected  for  study  in 


this  analysis  is  a  high  strength  steel  diaphragm  or  plate 
located  at  the  discrete  ring  at  frame  nine.  The  SAPIV 


finite  element  code  is  utilized  to  model  the  diaphragm.  Th 
diaphragm  has  been  divided  into  eight  "pie-shaped"  wedges 
which  are  further  subdivided  into  three  sections.  Twenty- 
four  plate  elements  are  used  in  this  model  with  twenty-five 
nodal  points.  Since  this  study  is  concerned  with  end-on 
loading,  the  node  at  the  center  of  the  plate  is  constrained 
to  move  only  in  the  x-direction  with  no  rotations.  The 
outer  nodes  (eighteen  through  twenty-five)  are  rigidly 
attached  to  the  shell.  This  is  accomplished  through  a 
modification  of  the  general  purpose  SAPIV  code  by  the 
developers  of  ELSHOK.  During  the  parametric  study,  the 
diaphragm  thickness  is  uniformly  varied  from  .25  inches  to 
9.25  inches.  Figure  5  is  a  TEKTRONIX  representation  of  the 
finite  element  model  of  the  substructure.  Appendix  B 
includes  one  of  the  SAPIV  input  codes  used  during  the 
parametric  study. 


IV.  ANALYSIS 


The  purpose  of  this  analysis  is  to  study  the  shell/ 
substructure  interaction  after  the  system  is  subjected  to  an 
end-on  load  while  varying  the  mass  and  stiffness  of  the 
internal  substructure.  Two  different  types  of  loading  are 
explored:  a  taper  charge  and  a  conventional  charge.  Taper 

charges  are  characterized  by  sustained  incident  pressure  on 
a  body  over  a  period  of  time.  They  are  used  in  underwater 
shock  tests  to  simulate  the  type  of  pressure  profile  which  a 
body  would  experience  from  a  nuclear  detonation.  Extensive 
use  of  scaling  is  involved  in  efforts  to  reproduce  the 
physical  behavior  of  an  object  in  response  to  a  nuclear 
charge  without  generating  the  extremely  high  pressures 
associated  with  actual  detonations.  Conventional  charges 
create  a  large  initial  pressure  pulse  followed  by  a  rapidly 
decaying  exponential.  ELSHOK  uses  inputs  of  charge  weight 
and  standoff  geometry  to  calculate  the  transient  velocities 
of  the  shell  and  substructure  in  response  to  these  loads. 


A.  ANALYSIS  PROCEDURE 

Once  the  shell  model  has  been  produced,  a  series  of 
SAPIV  runs  are  conducted  for  the  diaphragms  of  increasing 
thickness.  Each  run  generates  natural  frequencies  and 
accompanying  mode  shapes  for  each  model.  The  PICRUST 


V  V  -V 


program  reduces  the  data  obtained  through  SAPIV  into  a 
format  which  can  be  combined  with  the  shell-fluid  file  for 
processing  with  the  USLOB  code.  The  substructure  is 
connected  to  the  shell  during  this  phase  of  the  ELSHOK 
analysis.  The  diaphragm  is  attached  to  shell  segment  two  at 
frame  nine  (node  eleven) . 

The  weight,  type,  and  location  of  the  charge  is 
specified  in  the  USLOB  input  code.  To  set  up  an  end-on  load 
referenced  to  the  global  coordinate  system,'  the  charges  are 
placed  in  the  negative  X-direction  at  a  distance  of  70  feet 
(840  inches) .  This  distance  is  typical  of  stand-off 
distances  used  during  actual  tests.  For  the  taper  charge,  a 
typical  incident  pressure- time  history  shown  in  Figure  6  is 
used  to  simulate  the  shock  wave  loading.  In  order  to  make  a 
meaningful  comparison  between  the  structural  responses  from 
a  taper  charge  and  a  conventional  charge,  it  is  necessary 
that  they  have  equivalent  shock  wave  impulses.  In  the  field 
of  underwater  shock  analysis,  "impulse”  refers  to  the  time- 
integral  of  the  pressure  profile  [Ref.  9].  The  impulse  of 
unit  area  of  the  shock  wave  front  up  to  a  time  t  after  its 
arrival  is  given  by: 

I(t)  =  /  P(t)dt  (3) 

•'0 

The  impulse  of  the  taper  charge  is  found  by  calculating  the 
area  under  the  pressure-time  history  in  Figure  6.  As  was 


stated  earlier,  a  conventional  charge  follows  a  pressure¬ 
time  history  given  by  equation  (1): 

P(t)  «  Kj^(w^/^/R)  ^  exp(-t/6^)  (1) 

where : 

®o  “  (w^'^^/R)  ^  =  time  constant  of  exponential 

decay  (msec) 

HBX-1  is  a  popular  explosive  used  in  underwater  shock  tests, 
so  it  was  decided  to  use  it  in  this  study  for  the 
conventional  charge.  K2,  K^,  and  are  constants  which 

depend  on  the  type  of  explosive  used.  The  values  used  in 
this  study  are  [Ref.  11]: 

=  3.8354314E+05 

K2  =  1.14’4 

Kj  =  3.03131E-05 

K4  =  -.247 

By  matching  the  impulses  between  the  taper  charge  and  the 
conventional  charge,  a  charge  weight  of  352  lbs  of  HBX-1  is 
found  to  yield  an  equivalent  impulse.  The  peak  pressure 
generated  by  the  charge  is  1620  psi  which  creates  the 
incident  pressure-time  history  shown  in  Figure  7.  The 
conventional  charge  expends  most  of  its  impulse  in  the  first 


millisecond,  while  the  taper  charge  takes  over  four  mil 


seconds  to  expend  an  equivalent  amount  of  impulse. 

I 


Bec  = 


of  the  sharp  rise  in  pressure  experienced  from  a  conven 


195.44  Hz  for  the  conventional  charge.  A  series  of  SAPIV 
calculations  are  conducted  on  diaphragms  of  varying 
thickness  in  order  to  identify  the  fundamental  frequency  of 
each  substructure.  These  values  are  given  in  Table  II.  An 
analysis  is  performed  on  the  shell/substructure  system  for 
each  diaphragm  thickness  to  observe  the  coupled  velocity¬ 
time  history  response  as  the  substructure  mass  grows  in 
relation  to  that  of  the  shell. 


TABLE  II 


FUNDAMENTAL  FREQUENCIES  OF  DIAPHRAGM  SUBSTRUCTURE 


Thickness  (In) 


Frequency  (Hz) 


V.  RESULTS 


The  velocity  time-history  response  begins  at  time  zero 
when  the  incident  pressure  wave  from  the  underwater 
explosion  strikes  the  shell.  With  end-on  loading »  the 
forward  endplate  receives  the  initial  contact  with  velocity 
indicated  along  the  longitudinal  axis.  The  empty  shell 
response  is  analyzed  first,  followed  by  a  discussion  of  the 
coupled  effects  when  the  diaphragm  achieves  resonance. 

A.  VELOCITY-TIME  HISTORY  RESPONSE — EMPTY  SHELL 

The  empty  shell  velocity-time  history  responses  for  the 
conventional  and  taper  charges  are  contained  in  Figures  8 
through  20.  It  can  be  seen  that  the  conventional  charge 
causes  a  much  greater  transient  response  in  the  shell.  In 
both  cases,  the  velocity  is  rapidly  damped  by  interaction 
between  the  shell  and  the  infinite  fluid  medium.  It  is 
interesting  to  observe  the  differences  in  response  between 
the  various  locations  on  the  shell.  The  forward  endplate  is 
perpendicular  to  the  shock  wave,  receiving  the  full  force 
from  the  incident  pressure  pulse.  Inertial  forces  are 
generated  by  the  endplate  causing  severe  local  deformation 
in  the  nearest  nodes  on  the  cylinder  due  to  the  rigid 
coupling  between  segments.  The  motion  propagates  down  the 
shell  along  the  longitudinal  axis  until  it  reaches  the  aft 
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Velocity-Time  History  Response  (Empty  Shell) — Frame  Twenty-Three 


VELOCITY  RESPONSE  -  EMPTY  SHELL 

SEGMENT  TWO  -  NODE  THIRTY  SEVEN 
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Figure  17.  Velocity-Time  History  Response  (Empty  Shell) — Frame  Thirty-Six 


VELOCITY  RESPONSE  -  EMPTY  SHELL 

SEGMENT  TWO  -  NODE  FORTY  ONE 


Velocity-Time  History  Response  (Empty  Shell) — Frame  Forty 


endplate  over  a  millisecond  after  initial  shock  wave 
contact.  At  first  the  aft  endplate  moves  in  the  negative 
X-direction  as  the  "ripple"  proceeds  down  the  structure. 

Then  it  experiences  a  violent  transient  response  as  it  is 
affected  not  only  by  the  inertial  forces  from  the  shell,  but 
by  the  reflected  pressure  pulses  from  the  water  it 
displaces.  After  the  initial  shock,  the  shell  settles  into 
an  "accordion"  effect,  expanding  and  contracting  as  the 
forward  and  aft  endplates  move  in  opposite  directions. 

B.  VELOCITY-TIME  HISTORY  RESPONSE— COUPLED  SHELL/ 

SUBSTRUCTURE  SYSTEM 

From  the  results  of  the  analyses  performed  on  the  coupled 
shell/substructure  system.  Tables  III  and  IV  were  created  in 
order  to  identify  when  the  diaphragm  achieved  resonance.  Q/S 
is  the  frequency  ratio  between  the  fundamental  frequency  of 
the  simply-supported  diaphragm  (R)  and  the  dominant  frequency 
observed  in  the  velocity-time  history  response  of  the  sub¬ 
merged  empty  stiffened  shell  (R) .  M1/M2  is  the  mass  ratio 

between  the  substructure  (Mj^)  and  the  shell  with  endplates 
(M2) .  Figure  21  is  a  plot  of  the  velocity  of  the  diaphragm 
center  versus  the  frequency  ratio  R/S.  It  reveals  that  in 
both  the  taper  and  conventional  charge  responses,  the 
diaphragm  achieves  its  maximum  velocity  or  resonance  when  the 
fundamental  frequency  of  the  substructure  is  nearly  equal  to 


the  excitation  frequency  from  the  shell.  Resonance  occurs 
when  the  mass  ratio  is  approximately  .08.  Since 


TABLE  III 

MAXIMUM  VELOCITY  RESPONSE  OF  DIAPHRAGM  SUBSTRUCTURE 
WITH  VARYING  THICKNESS  (TAPER  CHARGE) 


Thickness  (In) 


M1/M2 


Velocity  (ft/sec) 


TABLE  IV 

MAXIMUM  VELOCITY  RESPONSE  OF  DIAPHRAGM  SUBSTRUCTURE 
WITH  VARYING  THICKNESS  (CONVENTIONAL  CHARGE) 


Thickness  (In) 


M1/M2 


Velocity  (ft/sec) 


Figure  22.  Two  Degrees  of  Freedom  System 


where:  Kl  *  overall  stiffness  of  the  substructure 
K2  »  overall  stiffness  of  the  shell 

it  can  be  said  that: 

Kl/Ml  =  K2/M2  =  K1/.08  M2 

or 

K1/K2  *  .08 

The  coupled  shell  and  substructure  can  then  be  modeled  as  a 
two  degrees  of  freedom  system  as  shown  in  Figure  22  where  Ml 
and  Kl  represent  the  substructure  and  M2  and  K2  represent 
the  shell.  The  large  mass  (M2)  has  an  uncoupled  natural 
frequency  (f2)  which  is  equal  to  <v/K2/M2’.  This  is 
equivalent  to  the  dominant  frequency  (ff)  of  the  empty  shell. 
Since  this  is  a  two  degrees  of  freedom  system,  it  will 
possess  two  natural  frequencies  for  the  coupled  system.  It 
has  been  found  in  coupled  systems  that  for  decreasing  mass 
ratio  (M1/M2)  and  for  increasing  stiffness  ratio  (K1/K2) , 
the  first  and  second  modal  frequencies  of  a  two  degrees  of 
freedom  system  approach  the  decoupled  frequencies  fj^  and  f2 
of  two  single  degrees  of  freedom  systems,  respectively  [Ref. 
12].  The  first  natural  frequency  is  bounded  by  f2»  and  fj^ 
is  bounded  by  the  second  natural  frequency.  From  Figure  23, 
it  is  seen  that  with  M1/M2  =  K1/K2  =  .08,  the  dominant 


frequency  of  the  shell  is  approximately  equal  to  the  first 
modal  frequency  of  the  system.  The  result  is  that  the 

shell/substructure  system  achieves  resonance  when  excited  at 
the  dominant  frequency  of  the  empty  shell. 

With  the  diaphragm  at  resonance.  Figures  24  and  25 

illustrate  that  the  transient  response  of  the  diaphragm  is 

much  more  pronounced  than  that  of  the  shell  due  to  the  large 

differences  between  the  two  masses  as  well  as  the  soft 

spring  between  them.  Initially  the  diaphragm  is  drawn  in 

the  negative  X-direction  by  inertial  forces  as  the  shell  is 

displaced  by  the  shock  wave.  With  the  increased  motion  of 

the  shell  at  frame  nine,  the  diaphragm  receives  more  and 

more  energy.  At  approximately  t  =  20  msec,  the  diaphragm  is 

at  its  maximum  velocity  while  the  shell  is  relatively  calm. 

At  this  point,  the  diaphragm  begins  to  transfer  energy  back 

2 

to  the  shell.  Since  kinetic  energy  is  equal  to  1/2  MV  ,  the 
diaphragm  is  a  small  mass  with  high  velocity  while  the  much 
larger  mass  of  the  shell  moves  at  a  lower  velocity.  The 
energy  of  the  system  is  transferred  between  the  shell  and 
substructure  as  the  overall  motion  is  damped  by  interaction 
between  the  shell  and  the  fluid  medium.  The  transfer  of 
energy  between  the  shell  and  the  diaphragm  creates  an  effect 
known  as  the  "beating"  phenomenon.  This  effect  is  seen  in 
the  velocity  time  history  responses  for  the  case  of 
diaphragm  resonance  in  Figures  26  through  39.  The  effect  is 


Velocity-Time  History  Response  (Resonant  Case) — Frame  Zero 


Velocity-Time  History  Response  (Resonant  Case) --Frame  Five 
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Velocity-Time  History  Response  (Resonant  Case) — Frame  Thirteen 
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igure  37.  Velocity-Time  History  Respi 


quite  pronounced  at  the  endplates,  because  they  act  as  "hard 
spots"  receiving  the  energy  transmitted  through  the  shell. 
The  cylindrical  portion  of  the  shell  acts  as  a  conduit  for 
the  energy  transfer.  Again,  the  nodes  on  the  shell  nearest 
the  endplate  experience  some  local  deformation  due  to  the 
rigid  coupling  between  segments. 

The  empty  shell  velocity-time  history  response  is 
compared  with  the  resonant  diaphragm  case  in  Figures  40 
through  45.  Of  particular  interest  is  the  fact  that  the 
velocity  profiles  are  virtually  identical  for  the  first 
fifteen  milliseconds.  During  this  period  of  time,  the  shell 
is  receiving  energy  from  the  shock  wave  and  transferring  a 
portion  of  it  to  the  diaphragm.  Because  of  the  differences 
in  mass  between  the  shell  and  substructure,  not  much  energy 
is  lost  from  the  shell.  The  peak  velocities  do  not  vary  by 
more  than  five  feet  per  second  during  this  time  period.  The 
difference  arises  at  t  =  35  msec  when  the  shell  receives 
energy  from  the  resonating  plate.  The  decaying  velocity 
response  is  disturbed  during  the  exchange  of  energy  which  is 
observed  in  the  beating  effect. 

As  the  thickness  of  the  diaphragm  is  increased  beyond 
the  resonant  condition,  the  interactions  between  the  shell 
and  substructure  are  reduced  dramatically.  The  shell/ 
substructure  velocity-time  history  remains  within  10%  of  the 
empty  shell  response  until  the  diaphragm  thickness  is  2.25 
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inches  thick  or  M1/M2  =  .164.  By  the  time  it  is  9.25  inches 
thick/  it  represents  over  half  the  mass  of  the  shell,  and  it 
acts  as  if  it  was  an  integral  part  of  the  shell  slaved  to 
the  motions  of  the  forward  endplate. 


VI.  CONCLUSIONS 


Upon  inspection  of  the  velocity-time  history  responses 
at  various  locations  along  the  shell,  one  can  see  why  the 
potential  for  damage  is  so  great  from  an  underwater 
explosion  off  the  bow  or  stern  of  a  submarine.  While  the 
cylindrical  portion  for  the  most  part  serves  to  transmit 
inertial  forces  with  little  velocity  response,  the 
extremities  of  the  hull  experience  violent  transient 
responses  which  may  be  passed  to  attached  internal  equipment 
located  in  those  areas.  The  inertial  influence  of  the  mass 
of  the  endplates  on  neighboring  nodes  of  the  shell  is 
apparent  in  both  the  empty  shell  and  resonant  cases.  The 
response  from  a  conventional  charge  is  seen  to  be  much  more 
pronounced  than  that  of  a  taper  charge  of  equivalent  impulse 
due  to  the  rapid  expenditure  of  energy  against  the  forward 
endplate . 

The  results  obtained  from  the  resonant  diaphragm  case 
support  the  theory  that  with  decreasing  mass  ratio  between 
substructure  and  shell,  the  first  and  second  modal 
frequencies  of  the  two  degrees  of  freedom  system  approach 
the  decoupled  frequencies  of  the  substructure  and  shell  in 
two  single  degree  of  freedom  systems.  This  allows  the 
motion  of  hull  to  drive  the  diaphragm  into  resonance  when 
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the  dominant  frequency  of  the  shell  motion  is  roughly  equal 
to  the  fundamental  frequency  of  the  substructure.  It  is 
possible  that  this  could  occur  during  an  underwater  shock 
test  of  a  small  piece  of  equipment  in  a  large  SSTV.  This 
might  result  in  the  failure  of  the  equipment  in  a  scenario 
which  would  not  occur  in  a  real  submarine  with  its  myriad  of 
internal  masses  and  substructures. 

The  large  difference  in  mass  between  the  substructure 
and  the  shell  results  in  a  tremendous  dynamic  amplification 
in  the  response  of  the  diaphragm  at  resonance.  As  energy  is 
transferred  between  the  two  bodies,  the  large  mass  of  the 
shell  causes  its  velocity  to  be  quite  small  compared  to  that 
of  the  diaphragm.  As  is  seen  in  Figures  40  through  45,  the 
initial  resonant  case  velocity  response  is  virtually 
unchanged  from  the  empty  shell  response.  If  the  mass  of  the 
internal  equipment  was  sufficiently  small  in  comparison  to 
an  SSTV  used  in  a  series  of  shock  tests,  one  ELSHOK  calcu¬ 
lation  could  be  performed  on  the  empty  shell  to  generate  a 
set  of  velocity  response  curves.  The  velocity-time  history 
response  at  the  attachment  point  could  then  be  input  into 
the  SAPIV  code  for  the  finite  element  model  of  any  number  of 
small  pieces  of  equipment,  and  the  substructure  transient 
response  could  be  obtained  without  having  to  run  a  separate 
ELSHOK  calculation  for  each  shell/substructure  system.  This 
would  result  in  considerable  savings  in  computer  time. 


Although  ELSHOK  is  a  powerful  tool  which  can  be  used  to 
accurately  predict  the  results  of  underwater  explosion  tests 
on  equipment,  some  shortfalls  exist  in  the  program  which 
limit  its  use  in  design  applications.  At  the  present  time, 
there  is  no  direct  method  to  take  the  effects  of  the  surface 
cut-off  wave  or  bottom  reflection  into  account.  The  surface 
cut-off  wave  forms  in  shallow  underwater  explosions  when  the 
compressive  shock  wave  reflects  off  of  the  surface  of  the 
water  causing  a  tensile  wave  to  propagate.  Bottom 
reflections  form  in  deep  explosions  as  the  shock  wave  is 
reflected  off  of  the  floor  of  the  ocean.  ELSHOK  assumes 
that  the  explosion  is  occurring  in  an  infinite  acoustic 
medium  with  a  single  shock  wave.  Another  weakness  in  the 
program  is  the  fact  that  it  cannot  process  oblique 
explosions;  charges  must  be  placed  in  a  pure  end-on  or 
side-on  configuration.  Shock  tests  conducted  in  accordance 
with  MIL-S-901D  are  either  end-on  or  side-on  loading, 
however,  so  this  is  not  a  serious  drawback  in  predicting  the 
transient  motion  of  equipment  in  SSTV's. 

It  is  recommended  that  further  study  be  conducted  in 
this  area  using  more  refined  shell  models  resembling  full- 
scale  submarines.  Once  a  shell  model  is  produced  for  a 
given  submarine  class,  researchers  can  create  different 
finite  element  models  of  actual  equipment  to  observe  the 
interactions  between  the  submarine  hull  and  equipment  of 


varying  mass  and  shape  when  subjected  to  end-on  or  side-on 
loads.  It  is  felt  that  ELSHOK  is  an  effective  method  to 
give  an  indication  of  the  response  of  equipment  in 
underwater  explosion  tests  without  needless  destruction  from 
unexpected  failures.  Its  continued  use  prior  to  UNDEX  tests 
will  lead  to  savings  in  both  time  and  effort. 
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APPENDIX  A 


PROGRAM  USED  TO  CONVERT  POSLOB  PUNCH  CARD  FILES 
TO  FORMAT  ACCEPTANCE  BY  EASYPLOT 

Along  with  the  TEKTRONIX  plots  of  the  velocity-time 
histories,  PUSLOB  generates  velocity  punch  card  files  for 
each  requested  node.  The  punch  card  file  consists  of  a 
block  of  velocities  in  exponential  format  the  size  of  which 
depends  upon  the  number  of  time  points  requested.  The 
EASYPLOT  program  developed  for  the  IBM  3033  by  Mr.  John 
Mainwaring  is  a  quick  interactive  method  to  produce  graphs 
using  DISSPLA  without  having  to  create  a  plotting  program. 

It  will  accept  up  to  100  data  points  in  tabular  form 
provided  it  is  in  a  column  format  for  X-Y  input.  The 
PLOTCONVl  code  was  created  to  accept  data  from  two  separate 
punch  card  files  for  taper  and  conventional  charge 
comparisons  and  combine  them  into  one  data  file  for  each 
node  for  input  into  EASYPLOT. 

In  order  to  prevent  data  overflow,  the  time  step  must  be 
sized  in  USLOB  and  PUSLOB  to  generate  no  more  than  IOC 
velocity  data  points.  Both  the  PLOTCONVl  VAX  11/780  VMS 
computer  command  file  and  the  FORTRAN  code  are  included  in 
this  appendix.  The  command  file  must  be  modified  prior  to 
each  run  to  specify  the  punch  card  files  being  read  and  the 
titles  of  the  new  data  files. 
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PLOTCONVl.COM  -  -  COMMAND  FILE  FOR  EXECUTION  OF  PLOTCONVl 


ASSIGN  PUSL48.PUN  FOROll 
ASSIGN  PUSL49.PUN  F0RQ12 
ASSIGN  MAXSHL1.DAT  FOR021 
ASSIGN  MAXSHL2.DAT  F0R022 
ASSIGN  MAXSHL3.DAT  F0R023 
ASSIGN  MAXSHL4.DAT  F0R024 
ASSIGN  MAXSHL5.DAT  F0R02S 
ASSIGN  MAXSHL6.DAT  FOR026 
ASSIGN  MAXSHL7.DAT  F0R027 
ASSIGN.MAXSHLa.DAT  FOR028 
ASSIGN  MAXSHL9.DAT  FOR029 
ASSIGN  MAXSHL10.DAT  F0R030 
ASSIGN  MAXSHL11.DAT  FQR031 
ASSIGN  MAXSHL12.DAT  F0R032 
ASSIGN  MAXSHL13.DAT  FQR033 
ASSIGN  MAXSHL14.DAT  F0R034 

RUN  PLOTCONVl 


C 

C 

C 

C 

C 


PROGRAM  PLOTCONVl 

program  IS  TO  READ  TWO  PUNCH-CARD 
BE  ACCEPTED  BY  EASYPLOT®^''^'^^  ^  WHICH  CAN 

CHARACTERXIO  LABEL 

CHARACTERX40  LTITLE 

COMMON/CSOLVE/DTRECS 

NCURV=14 

DO  10  K»1,NCURV 

5|AD( 11 , 999 )  LTITLE( K) , NRECS. DTRECS 
READC11,998)  LABELCK),  TAGdC) 

READai,997)  (VELCJ.K),  J  =  l, NRECS) 

R|AD(12.999)  LTITLE(K). NRECS, DTRECS 
READ(12,998)  LABEL(K) ,TAG(IO 
1,  ""“'■'•’"r  jn.NRECS) 

999  F0RMAT(A40,I5,1PE11.4) 

998  FORMAT(A10,F10.5) 

997  F0RMAT(1P7E11.4) 

DO  20  L=1,NCURV 

T=0.0 

LL=L+20 

DO  30  M=l, NRECS 
VEL0<M,1)=T 
DTRECS1=1000.0XDTRECS 
T=T+DTRECS1 

VEL0(M,2)=VEL(M,L)/12.0 
VELO(M,3)=VEL1(M,L)/12.0 
30  CONTINUE 

DO  40  KK=1, NRECS 

40  CO^TINuk'”'’  (VEL0<KX,JJ),JJ=1,3) 

W5ITE(LL,995)  LTITLEIL), NRECS, DTRECS 
WRITE(LL,994)  LABEL(L) ,TAG<L) 

20  CONTINUE 
996  F0RMAT(3F1S.5) 

995  FORMAT(A40,I5,1PE11.4) 

994  FORMAT(A10,F10.S) 

STOP 
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APPENDIX  B 


TYPICAL  ELSHOK  INPUT  CODES — SHELL  AND  SUBSTRUCTURE 

The  input  codes  which  follow  are  for  the  resonant 
diaphragm  case  for  both  taper  and  conventional  charges.  Th 
only  differences  between  the  two  is  in  the  USLOB  code  where 
the  type  of  charge  is  specified.  All  the  input  codes  are 
explained  in  detail  in  the  ELSHOK  users  manual  [Ref.  4]. 

1.  B0S0R4  INPUT  DA'^A 

The  B0S0R4  input  code  is  created. by  working  inter¬ 
actively  with  the  computer.  The  prompt  following  the  $ 
symbol  appears,  and  the  user  supplies  the  required 
information.  The  three  segments  (two  endplates  and 
cylinder)  are  modeled  separately,  and  then  they  are 
"connected"  in  the  global  data  section.  A  separate  input 
code  is  required  for  each  circumferential  harmonic  included 
in  the  calculation.  B0S0R4  output  provides  the  user  with 
the  in-vacuo  free-fr.ee  modes  and  natural  frequencies  of  the 
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INDIC  «  ANALTSIS  TYPE  INDICATOR 

NPRT  a  OUTPUT  OPTIONS  (laMIHINUM.  2aNE0IUM.  SaHAXINUM) 

ISTRESa  OUTPUT  CONTROL  (OaRESULTANTS.  laSXGMA,  2aEPSIL0N) 

NSEO  a  number  OF  SHELL  SEGMENTS  (LESS  THAN  23) 

SEGMENT  NUMBER  11111111 
NMESH  a  number  OP  NODE  POINTS  (3  ■  MIN. t  98  a  nAX. >(  1) 

NTYPEHa  CONTROL  INTEGER  (1  OR  2  OR  3)  FOR  NODAL  POINT  SPACING 
NSHAPEa  INDICATOR  (1.2  OR  A)  FOR  GEOMETRY  OF  MERIDIAN 
R1  a  radius  at  beginning  of  segment  (SEE  P.  66) 

21  a  axial  coordinate  AT  BEGINNING  OF  SEGMENT 
R2  a  RADIUS  AT  END  OF  SEGMENT 

22  a  axial  coordinate  at  end  of  segment 

IMP  a  indicator  for  imperfection  (OaNONE,  l-SOME) 

NTYPEEa  CONTROL  (1  OR  3)  FOR  REFERENCE  SURFACE  LOCATION 
2VAL  a  distance  FROM  LEFTMOST  SURF.  TO  REFERENCE  SURF. 

DO  YOU  NANT  TO  PRINT  OUT  R(S).  R*(S>.  ETC.  FOR  THIS  SEOMENTT 
NRINOSa  number  (MAXaEO)  OF  DISCRETE  RINGS  IN  THIS  SEGMENT 
RaELASTIC  FOUNDATION  MODULUS  (E.G.  LB/INNNS>IN  THIS  SEO. 
LINTYPa  INDICATOR  (0.  1.  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 
NLTYPEaCONTROL  (0.1. 2.3)  FOR  TYPE  OF  SURFACE  LOADING 
NHALL«INDEX  (1.  2.  A.  3.  6.  7.  8)  FOR  HALL  CONSTRUCTION 
E  a  YOUNG'S  MODULUS  FOR  SKIN 

U  a  POISSON'S  RATIO  FOR  SKIN 

SM  aHASS  DENSITY  OF  SKIN  (E.G.  ALUM. a. 00023  LB-SEC3i«2/IHmA) 
ALPHA  a  COEFFICIENT  OF  THERMAL  EXPANSION 

ANRS  a  CONTROL  (0  OR  1)  FOR  ADDITION  OF  SHEARED  STIFFENERS 
SUR  a  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -1) 

DO  YOU  HANT  TO  PRINT  OUT  THE  C(I.J)  AT  MERIDIONAL  STATIONS? 

DO  YOU  HANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIANT 

SEGMENT  NUMBER  22222222 
NMESH  a  NUMBER  OF  NODE  POINTS  (3  a  MIN.i  98  ■  NAX.)(  2) 

NTYPEHa  CONTROL  INTEGER  (1  OR  2  OR  3)  FOR  NODAL  POINT  SPACING 
NSHAPEa  INDICATOR  (1.2  OR  A)  FOR  GEOMETRY  OF  MERIDIAN 
R1  a  radius  at  beginning  OF  SEGMENT  (SEE  P.  66) 

21  a  AXIAL  COORDINATE  AT  BEGINNING  OF  SEGMENT 
R2  a  radius  at  end  OF  SEGMENT 

22  a  axial  coordinate  AT  END  OF  SEGMENT 

IMP  a  indicator  for  imperfection  (OaNONE.  laSOME) 

NTYPEZa  CONTROL  (1  OR  3)  FOR  REFERENCE  SURFACE  LOCATION 
2VAL  a  distance  FROM  LEFTMOST  SURF.  TO  REFERENCE  SURF. 

DO  YOU  HANT  TO  PRINT  OUT  R(S).  R'(S).  ETC.  FOR  THIS  SEOMEHTt 
NRINGS*  NUMBER  <NAXa20)  OF  DISCRETE  RINGS  IN  THIS  SEGMENT 

CONTROL  FOR  IDENTIFICATION  OF  RING  LOCATION  (2a2.  3aR) 
AXIAL  COORDINATE  OF  ITH  RING,  2(  1) 

AXIAL  COORDINATE  OF  ITH  RING. 

AXIAL  COORDINATE  OF  ITH  RING. 

AXIAL  COORDINATE  OF  ITH  RING. 

AXIAL  COORDINATE  OF  ITH  RING. 

AXIAL  COORDINATE  OF  ITH  RING, 


NTYPE 
2(1) 

2(1) 

2(1) 

2(1) 

2(1) 

2(1) 

HTYPERa  TYPE  (0  OR 
NTYPERa  type  (0  OR 
NTYPERa  TYPE  (0  OR 
HTYPERa  TYPE  (0  OR 
NTYPERa  TYPE  (0  OR 
NTYPERa  TYPE  (0  OR 


2(  2) 

2(  3) 

2(  A) 

2(  5) 

2(  6) 

1  OR  2  OR  A  OR  5)  OF  DISCRETE  RING  NO.(  1) 

1  OR  2  OR  A  OR  5)  OF  DISCRETE  RING  NO.(  2) 

1  OR  2  OR  A  OR  5)  OF  DISCRETE  RING  NO.(  3) 

1  OR  2  OR  A  OR  5)  OF  DISCRETE  RING  NO.(  A) 

1  OR  2  OR  A  OR  5)  OF  DISCRETE  RING  NO.(  5) 

1  OR  2  OR  A  OR  5)  OF  DISCRETE  RING  NO.(  6) 


E  a  YOUNG'S  MODULUS  OF  RINO(  1) 

A  a  CROSS  SECTION  AREA  OF  RINO(  1) 

lY  a  moment  of  INERTIA  ABOUT  Y-AXIS  (SEE  FIG.  ON  P.70X  1) 

IX  a  moment  of  inertia  ABOUT  X-AXIS(  1) 

IXY  a  PRODUCT  OF  INERTIAL  1) 

El  a  radial  component  of  ring  ECCENTRICITY  (SEE  P.  70)(  1) 

E2  a  axial  component  OF  RING  ECCENTR1CITY(  1) 

OJ  a  TORSIONAL  RIOtOITY(  1) 

RM  a  RIHO  material  DENSITY  (E.G.  ALUMINUMa . 0002333) (  1) 

E  a  YOUNG'S  MODULUS  OF  R1NC(  2) 

A  a  CROSS  section  AREA  OF  RING!  2) 

lY  a  moment  of  inertia  ABOUT  Y-AXIS  (SEE  FIG.  ON  P.70)(  2) 

IX  a  moment  of  inertia  about  X-AX1S(  2) 

IXY  a  PRODUCT  OF  INERTIAL  2) 
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FILE:  lOSSO 


DATA 


A1 


-Z.OtZSOO 
O.OOOOOOOE'^00 
0.9249000E-»'Ot 
0.7330000E>«3 
0.3000000E-»'Ot 
'  S.2SQ000 
11.&9S00 
2.750000 
O.OOOOOOOE-VOO 
-2.0A2SOO 
O.OOOOOOOEt^OO 
0.92A9000E-»OS 
0.7330000E>03 
0.3000000E'»OS 
0. 250000 
11.09000 
2.750000 
o.oooooooE>oo 
-2.002500 
O.OOOOOOOE^OO 
0.9209000E>0S 
0.73300006-03 
0.3000000E>aS 
0.250000 
11.09000 
2.750000 
O.OOOOOOOE'fOO 
-2.002500 
a.aooooooE>oo 
0.9249a00E4'OO 
0.73300006-03 
0.30000006>00 
7.125000 
7.535500 
2.375000 
0.00000006000 
-1.701250 
0.00000006000 
0.75016000 
0.73300006-03 
0.00000006000 
0 
0 
2 

0.30000006000 

0.3000000 

0.73300006-03 

0.0 

1 

1 

N 

Y 

0 

0.30000006000 

0.3000000 

0.73300006-03 

Y 

Y 

5.625000 

0.2500000 

1.250000 

N 

N 

H 

H 

10 

3 

1 

16.01250 

0.0 

o.aooaooo6«oo 

0.0 


*61  ■  RADIAL  COMFONENT  OF  RING  ECCENTRICITY  (SEE  F.  70)(  2) 

•62  ■  AXIAL  COMFONENT  OF  RINO  ECCENTRICITY!  2) 

•  OJ  >  TORSIONAL  RIGIDITY!  2) 

•  RM  •  RINO  NATERIAL  DENSITY  !E.O.  ALUMINUM-. 00025351!  2) 

«  E  •  YOUNG'S  MODULUS  OF  RINO!  3) 

•A  ■  CROSS  SECTION  AREA  OF  RINO!  3) 

•  lY  «  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  !SEE  FIO.  ON  F.70)!  3) 

•IX  ■  MOMENT  OF  INERTIA  ABOUT  X-AXIS!  3) 

•  IXY  ■  PRODUCT  OF  INERTIA!  31 

•61  ■  RADIAL  COMFONENT  OF  RINO  ECCENTRICITY  !SEE  F.  701!  31 

•62  •  AXIAL  COMFONENT  OF  RINO  ECCENTRICITY!  31 

•  OJ  >  TORSIONAL  RIGIDITY!  31 

•  RM  «  RING  MATERIAL  DENSITY  !E.O.  ALUMINUM-. 00025351!  31 

•  E  -  YOUNG'S  MODULUS  OF  RINO!  6) 

•A  -  CROSS  SECTION  AREA  OF  RINO!  61 

•  lY  -  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  !SEE  FIO.  ON  P.701!  61 

•IX  -  MOMENT  OF  INERTIA  ABOUT  X-AXIS!  61 

•IXY  -  PRODUCT  OF  INERTIA!  61 

•El  -  RADIAL  COMPONENT  OF  RINO  ECCENTRICITY  !SEE  P.  701!  61 

•  E2  -  AXIAL  COMPONENT  OF  RINO  ECCENTRICITY!  61 

•  OJ  -  TORSIONAL  RIGIDITY!  61 

«  RM  -  RINO  MATERIAL  DENSITY  lE.O.  ALUMINUM-. 00025351 I  61 

•  E  -  YOUNO'S  MODULUS  OF  RIM!  51 

•  A  «  CROSS  SECTION  AREA  OF  RING!  51 

•  lY  -  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  !SEE  FIO.  ON  P.701!  51 

•  IX  -  MOMENT  OF  INERTIA  ABOUT  X-AXIS!  5) 

•  IXY  -  PRODUCT  OF  INERTIA!  51 

•  El  -  RADIAL  COMPONENT  OF  RING  ECCENTRICITY  !SEE  P.  701!  51 

•  E2  >  AXIAL  COMPONENT  OF  RINO  ECCENTRICITY!  51 

•  OJ  -  TORSIONAL  RIGIDITY!  51 

•  RM  «  RINO  MATERIAL  DENSITY  !E.O.  ALUMINUM-. 00025351!  51 

•  E  -  YOUNG'S  MODULUS  OF  RING!  61 

•  A  -  CROSS  SECTION  AREA  OF  RING!  61 

•  lY  -  MOMENT  OF  INERTIA  ABOUT  Y-AXIS  !SEE  FIO.  ON  P.701!  61 

•  IX  •  MOMENT  OF  INERTIA  ABOUT  X-AXIS!  61 

•  IXY  «  PRODUCT  OF  INERTIA!  61 

•  El  •  RADIAL  COMPONENT  OF  RINO  ECCENTRICITY  !SEE  P.  701!  61 

•  E2  •  AXIAL  COMPONENT  OF  RINO  ECCENTRICITY!  61 

•  OJ  •  TORSIONAL  RIGIDITY!  61 

•  RM  •  RING  MATERIAL  OtOtStrY  tg.O.  ALUMINUM-. 00025351!  61 

•  K-EIASTIC  FOUNDATION  MODULUS  !E.O.  LB/INK»3>IN  THIS  SEO. 

•  LINTYP*  INDICATOR  !0.  I,  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 

•  NLTYPE-CONTROL  <0,1.2.31  FOR  TYPE  OF  SURFACE  LOADING 

•  NHALL-INOEX  (1.  2.  6.  5.  6.  7,  B1  FOR  HALL  CONSTRUCTION 

•  E  •  YOUNG'S  MODULUS  FOR  SKIN 

•  U  «  POISSON'S  RATIO  FOR  SKIN 

•  SM  -NAS*  DENSITY  OF  SKIN  !E.O.  ALUM. -.00025  IB-SECXN2/INKN61 

•  ALPHA  -  COEFFICIENT  OF  THERMAL  EXPANSION 

•  ANRS  -  CONTROL  <0  OR  11  FOR  ADDITION  OF  SMEARED  STIFFENERS 

•  SUR  -  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -11 

ARE  THERE  STRINGERS  !PLEASE  ANSMER  Y  OR  NIT 
ARE  THERE  RINGS  <PLEASE  ANSHER  Y  OR  NIT 
K2  -CONTROL  <0  OR  1}  FOR  INTERNAL  OR  EXTERNAL  RINGS 
E2  -  RINO  MODULUS 

U2  «  RINO  POISSON  RATIO 

ROMO-  RINO  NASS  DENSITY 

IS  THE  RINO  CROSS  SECTION  CONSTANT  IN  THIS  SEGMENTT 
IS  THE  RINO  CROSS  SECTION  RECTANGULAR  (Y  OR  NIT 
D2  -  ARC  LENGTH  BETWEEN  ADJACENT  RINGS  CCONSTANTl 
T2  -  THICKNESS  OF  RINO  !CONSTANT> 

H2  -  HEIGHT  OF  RINO  (CONSTANT) 

DO  YOU  HANT  TO  PRINT  OUT  THE  Cd.Jl  AT  MERIDIONAL  STATIONS! 
DO  YOU  WANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIAN! 


SEGMENT  NUMBER  33333333 
NMESH  -  NUMBER  OF  NODE  POINTS  (5  «  MIN.;  9S  -  MAX.)!  3) 
NTYPEH-  CONTROL  INTEGER  (1  OR  2  OR  3)  FOR  NODAL  POINT  SPACING 
NSHAPE-  INDICATOR  (1.2  OR  6)  FOR  GEOMETRY  OF  MERIDIAN 
R1  -  RADIUS  AT  BEGINNING  OF  SEGMENT  (SEE  P.  66) 

21  -  AXIAL  COORDINATE  AT  BEGINNING  OF  SEGMENT 
R2  -  RADIUS  AT  END  OF  SEGMENT 

22  -  AXIAL  COORDINATE  AT  END  OF  SEGMENT 
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FILEt  lOSSO 


DATA 


A1 


0 

3 

2.500000 

N 


o.iooooooe-»o< 

0.3200000 

0.253SOOOE-03 

0.0 


-I 

0 

0 

0 

1 

30 

0 . 0OOOOOOE-»OO 
o.ooQaaoaE>oo 
O.OOOOOOOE-^OO 
1 


1 

1 

1 

10 

1 

1 

1 

1 

o.ooooooaE>oo 

o.oooaoooE>oo 

Y 

1 


10 

0 


1 

1 

2 

AS 

1 

1 

1 

1 

O.OOOOOOOE+OO 

O.OOOOOOOE+00 

Y 


IMP  «  INDICATOR  FOR  INFERFECTIOH  (O-NONE.  1«S0ME) 

NTYFEZ>  CONTROL  (1  OR  3)  FOR  REFERENCE  SURFACE  LOCATION 
2VAL  ■  DISTANCE  from  LEFTMOST  SURF.  TO  REFERENCE  SURF. 

DO  YOU  HANT  TO  PRINT  OUT  RCS).  R*(S).  ETC.  FOR  THIS  SEGMENT? 
NRINOSa  NUMBER  (MAX>20)  OF  DISCRETE  RINGS  IN  THIS  SEGMENT 
K«ELASTIC  FOUNDATION  MODULUS  (E.O.  LB/IHn3)IN  THIS  SEG. 
LINTYP«  INDICATOR  (0.  1.  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 
NLTYPE«CONTROL  (0.1.2.31  FOR  TYPE  OF  SURFACE  LOADING 
NHALL-INDEX  (1.  2.  4.  S.  «.  7.  0)  FOR  HALL  CONSTRUCTION 
E  ■  YOUNG'S  MODULUS  FOR  SKIN 

U  ■  POISSON'S  RATIO  FOR  SKIN 

SM  -NASS  DENSITY  OF  SKIN  (E.O.  ALUM. -.00025  LB-SECXH2/INNH4} 
ALPHA  -  COEFFICIENT  OF  THERMAL  EXPANSION 

ANRS  -  CONTROL  (0  OR  1)  FOR  ADDITION  OF  SMEARED  STIFFENERS 
SUR  -  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -1) 

DO  YOU  HANT  TO  PRINT  OUT  THE  C(I.J)  AT  MERIDIONAL  STATIONS? 
DO  YOU  HANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIAN? 


GLOBAL  DATA.... 

NLAST  «  PLOT  OPTIONS  ('l-NONE.  t-GEONETRY.  1-U.V.H) 

STARTING  NUMBER  OF  CIRC.  HAVES  (BUCKLING  ANALYSIS) 
MINIMUM  NUMBER  OF  CIRC.  HAVES  (BUCKLING  ANALYSIS) 
MAXIMUM  NUMBER  OF  CIRC.  HAVES  (BUCKLING  ANALYSIS) 
INCREMENT  IN  NUMBER  OF  CIRC.  HAVES  (BUCKLING) 
NUMBER  OF  EIGENVALUES  FOR  EACH  HAVE  NUMBER 
PRESSURE  OR  SURFACE  TRACTION  MULTIPLIER 
I  TEMPERATURE  RISE  MULTIPLIER 
ANGULAR  VEL.  ABOUT  AXIS  OF  REVOLUTION  (RAD/SEC) 


NOB 

NMINB 

NMAXB 

INCRB 

NVEC 

P 

TEMP 

OMEGA 


NUMBER  OF  POLES  (PLACES  HHERE  R-0)  IN  SEGMENT!  1) 


CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  1  1  1  1 

IPOLE  -  NODAL  POINT  NUMBER  OF  POLE.  IPOLE(  1) 

AT  HOH  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 
IS  THIS  SEGMENT  JOINED  TO  ANY  LONER-NUMBERED  SEGMENTS? 
NUMBER  OF  POLES  (PLACES  HHERE  R-0)  IN  SEGMENT!  2) 


CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  2  2  2  2 

AT  HOH  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 

IS  THIS  SEGMEHT  JOINED  TO  ANY  LOHER-NUMBEREO  SEGMENTS? 

AT  HOH  MAY  STATIONS  IS  THIS  SEGMENT  JOINED  TO  PREVIOUS  SEGS.? 
INODE  «  NODE  IN  CURRENT  SEGMENT  (ISEO)  OF  JUNCTION.  INODE!  1) 
JSEO  «  SEGMENT  NO.  OF  PREVIOUS  SEGMENT  INVOLVED  IN  JUNCTION 
JNOOE  -  NODE  IN  PREVIOUS  SEONNT  (JSEG)  OF  JUNCTION 
lUSTAR-  AXIAL  DISPLACEMENT  (O-NOT  SLAVED.  1-SLAVED) 

IVSTAR-  CIRCUMFERENTIAL  DISPLACEMENT  (O-NOT  SLAVED.  1-SLAVED) 
IHSTAR-  RADIAL  DISPLACEMENT  (D-NOT  SLAVED.  1-SLAVED) 

ICHI  «  MERIDIONAL  ROTATION  (O-NOT  SLAVED.  1-SLAVED) 

D1  -  RADIAL  COMPONENT  OF  JUNCTURE  GAP 
02  -  AXIAL  COMPONENT  OF  JUNCTURE  CAP 

IS  THIS  CONSTRAINT  THE  SANE  FOR  BOTH  PREBUCKLINO  AND  BUCKLING? 
NUMBER  OF  POLES  (PLACES  HHERE  R-O)  IN  SEGMENT!  3) 


CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  3  3  3  3 

IPOLE  >  NODAL  POINT  NUMBER  OF  POLE.  IPOLE(  1) 

AT  HOH  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 

IS  THIS  SEGMENT  JOINED  TO  ANY  lOHER-NUMBERED  SEGMENTS? 

AT  HOH  NAY  STATIONS  IS  THIS  SEGMENT  JOINED  TO  PREVIOUS  SEGS.? 
INOOE  >  NODE  IN  CURRENT  SEGMENT  (I5EG)  OF  JUNCTION.  INODE!  1) 
JSEO  -  SEGMENT  NO.  OF  PREVIOUS  SEGMENT  INVOLVED  IN  JUNCTION 
JNODE  *  NODE  IN  PREVIOUS  SEGMNT  (JSEO)  OF  JUNCTION 
lUSTAR-  AXIAL  OISPIACENENT  (O-NOT  SLAVED.  1 -SLAVED) 

IVSTAR-  CIRCUMFERENTIAL  DISPLACEMENT  (O-NOT  SLAVED.  1-SLAVED) 
IHSTAR-  RADIAL  DISPLACEMENT  (O-NOT  SLAVED.  1-SLAVED) 

ICHI  -  MERIDIONAL  ROTATION  (O-NOT  SLAVED.  1-SLAVED) 

D1  «  RADIAL  COMPONENT  OF  JUNCTURE  GAP 
D2  *  AXIAL  COMPONENT  OF  JUNCTURE  GAP 

IS  THIS  CONSTRAINT  THE  SAME  FOR  BOTH  PREBUCKLINO  AND  BUCKLING? 
DO  YOU  HANT  TO  LIST  OUTPUT  FOR  SEGMENT!  1) 

00  YOU  HANT  TO  LIST  OUTPUT  FOR  SEGMENT!  2) 

DO  YOU  HANT  TO  LIST  OUTPUT  FOR  SEGMENT!  3) 

DO  YOU  HANT  TO  LIST  PREBUCKLINO  RESULTANTS  AND  RING  FORCES? 
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2.  ACESNID  INPUT  DATA 


The  ACESNID  computer  code  determines  the  virtual  mass 
array,  a  quantity  producing  the  late-time  contribution  of 
the  DAA.  Only  one  ACESNID  calculation  is  required  for  all 
circumferential  harmonics  of  the  shell.  The  bands  dividing 
the  areas  of  the  shell  are  identified  in  the  NLEFT,  NCYL, 
and  NRITE  entries,  while  the  fluid  properties  are  specified 
in  the  RHOFL  and  VSOUND  row. 


FILE>  ACEl  DATA  A1 


VIRTUAL  MASS  FOR  RING-STIFFENED  CYLINDER  HITH  FLAT  ENDS,  N«0,l 

mstart,nfinis,nfreq,nvmass.ncheck 

15  51  IS  30/  NLEFT.NCYL, NRITE. NHBOSG.NSYMF 

i  ®  ^  morder.nfends.nfcent.nfcmpt.nomit 

9.59684E-5  5,833E4  0,005  /  RHOFL, VSOUND, ERR 
1111110111/  OUTPUT  FLAGS 
CPS<1) 

12  /  sefs,jsgbeg,jptbeo,jsgend,jptend,left 

SEFS, JSGBEG, JPTBEG, JSGEND, JPTEND, CENTRAL 
SEFS, JSGBEG, JPTBEG, JSGEND, JPTEND, RIGHT 
SOURCES, JSGBEG, JPTBEG, JSGEND, JPTEND, LEFT 
SOURCES, JSGBEG, JPTBEG, JSGEND, JPTEND, CENTRAL 


1.0 
1  1 


47 

12 

12 

47 


3.  PIFLASH  INPOT  DATA 


The  PIFLASH  code  takes  data  from  the  shell  mode  files 
from  BOSOR4  and  the  virtual  mass  file  from  ACESNID  and 
combines  them  to  create  a  "shell-fluid  file"  containing  all 
the  information  concerning  the  shell  and  its  acoustic  medium 
for  further  calculations.  Since  only  the  breathing  modes 
are  activated  in  end-on  loading,  all  torsional  modes  are 
dropped  from  the  shell  file.  This  is  accomplished  in  the 
NJUSE  rows  where  the  modes  of  Table  I  are  identified. 


FILE:  PIEl  DATA  A1 


20000/  NUMBER, MTORSM.MPTW.MSYMS.MSYMP 
5  3  3  3  /  <MWETSG<K),NUSESG(IC),Ksl,2) 

336.4  /  GRAVITY 

SHELLOO  /  SMF,  N»0 

SHELLOl  /  SMF,  N»1 

19  2  /  <NJUSE<J),J»1,MITEMS) 

123/  (KORSG(K),K=1,NKORSG) 

2  4  6  7  9  12  13  14  19  20  21  22  23  25  26  27  28  29  30  /  JUSE,  N=0 
12/  JUSE,  N=l 


4.  SAPIV  INPUT  DATA 

The  substructure  is  modeled  using  the  standard  SAPIV 
code  with  slight  modifications  by  the  developers  of  ELSHOK. 
The  diaphragm  is  modeled  using  twenty-four  plates  and 
twenty-five  nodes.  The  location  of  each  node  is  identified 
in  cartesian  coordinates.  Each  plate  is  1.06  in  thick  for 
the  resonant  case.  Nodes  eighteen  through  twenty-five  are 
rigidly  attached  to  the  shell  in  the  last  eight  lines  of  the 
code.  The  SAPIV  code  determines  the  fixed-base  modes  and 
corresponding  natural  frequencies  as  well  as  the  uncon¬ 
strained  mass  and  stiffness  matrices  for  the  substructure. 


FILE:  SAP23  DATA  A1 


IDEALIZED  INTERNAL  EQUIPMENT  AT  FRAME  9 
25  10  10  100000/  NUMNP.NELTYP, . . .,N10SV,NRIGID 
0  1  0  1  1  1  1  1  0.0  0.0  0.0  0  0.0 
0  2  0  0  0  1  0  0  0.0  0.0  5.604  0  0.0 
03000100  0.0  3.963  3.963  0  0.0 
04000100  0.0  5.604  0.0  0  0.0 
05000100  0.0  3.963  >3.963  0  0.0 
0  6  0  0  0  1  0  0  0.0  0.0  -5.604  0  0.0 
07000100  0.0  -3.963  -3.963  0  0.0 
08000100  0.0  -5.604  0.0  0  0.0 
09000100  0.0  -3.963  3.963  0  0.0 
0  10  0  0  0  1  0  0  0.0  0.0  11.208  0  0.0 

0  11  0  0  0  1  0  0  0.0  7.925  7.925  0  0.0 

0  12  0  0  0  1  0  0  0.0  11.208  0.0  0  0.0 

0  13  0  0  0  1  0  0  0.0  7.925  -7.925  0  0.0 

0  14  0  0  0  1  0  0  0.0  0.0  -11.208  0  0.0 

0  15  0  0  0  1  0  0  0.0  -7.925  -7.925  0  0.0 

0  16  0  0  0  1  0  0  0.0  -11.208  0.0  0  0.0 

0  17  0  0  0  1  0  0  0.0  -7.925  7.925  0  0.0 

0  18  0  0  0  1  0  0  0.0  0.0  16.813  0  0.0 

0  19  0  0  0  1  0  0  0.0  11.888  11.888  0  0.0 

0  20  0  0  0  1  0  0  0.0  16.813  0.0  0  0.0 

0  21  0  0  0  1  0  0  0.0  11.888  -11.888  0  0.0 

0  22  0001000. 00. 0  -16.813  0  0.0 
0  23  0  0  0  1  0  0  0.0  -11.888  -11.888  0  0.0 

0  24  0  0  0  1  0  0  0.0  -16.813  0.0  0  0.0 

0  25  0  0  0  1  0  0  0.0  -11.888  11.888  0  0.0 

6  24  1  /  TYPE  6  (PLATE  AND  SHELL  ELEMENTS) 

1  .733E-03  0.0  0.0  0.0 

32.967E+06  9.89E+06  0.0  32.967E+06  0.0  11.538E+06 
0.0  0.0  0.0  0.0 
0.0  0.0  0.0  0.0 
0.0  0.0  0.0  0.0 
0.0  0.0  0.0  0.0 

0.0  0.0  0.0  0.0  /  ELEMENT  LOAD  MULTIPLIERS 
118900101. 06  000 
219200101. 06  000 
312300101. 06  000 
413400101. 06  000 
514500101. 06  000 
615600101. 06  000 
716700101. 06  000 
817800101. 06  000 
9  2  10  11  3  0  1  1  1.06  000 

15  8  16  17  9  0  1  0  1.06  000 

16  9  17  10  2  0  1  0  1.06  000 

17  10  18  19  11  0  1  1  1.06  000 

23  16  24  25  17  0  1  0  1.06  0  0  0 

24  17  25  18  10  0  1  0  1.06  0  0  0 

0  0  0.0  0.0  0.0  0.0  0.0  0.0  /  CONC  MASSES 
0.0  0.0  0.0  0.0  /  ELEMENT  LOAD  MULTIPLIERS 
0000. 00. 000. 0/  IFPR,IFSS, .. .,NFO, SHIFT 
80/  NFACE,NSYMTY 

18  1  1  1  1  0  0 

19  1  1  1  1  0  0 

20  1  1  1  1  0  0 

21  1  1  1  1  0  0 

22  1  1  1  1  0  0 

23  1  1  1  1  0  0 

24  1  1  1  1  0  0 

25  1  1  1  1  0  0 

END  OF  DATA  (DUMMY  TITLE) 


5.  PICROST  INPUT  DATA 


The  PICRUST  code  takes  data  from  the  substructure  mode 
file  from  SAPIV  and  reorganizes  it  to  facilitate  the 
solution  of  the  transient  response  equations.  The  location 
of  the  substructure  attachment  points  to  the  shell  is 
identified  for  each  connecting  node. 


FILE!  PICRUSTl  DATA  A1 


1111011100 

0111110101/  OUTPUT  FLAGS  1-20 
200/  NJUSE.NHNSOB.NHMBAR 
18  2  11  0.0 

19  2  11  45.0 

20  2  11  90.0 

21  2  11  135.0 

22  2  11  180.0 

23  2  11  225.0 

24  2  11  270.0 

25  2  11  315.0  /  NIPSUB,LBOSEG,LBOSPT,AMGDEG 
0.0  /  DEGROT 


6.  OSLOB  INPUT  DATA 

The  two  USLOB  input  codes  for  the  taper  and  conventional 
charges  appear  on  the  following  page.  USL46  DATA  is  the 
code  for  the  taper  charge.  It  contains  nine  discrete  data 
points  from  the  pressure-time  history  in  seconds  and  pounds 
per  square  inch.  There  are  1600  time  points  used  in  the 
numerical  integration  of  the  governing  equations  with 
information  saved  every  sixteenth  point,  so  only  100 
velocity-time  data  points  are  retained  for  EASYPLOT  display. 
The  charge  is  placed  along  the  negative  X-axis  at  840 
inches.  USL47  DATA  contains  the  information  for  the 
conventional  charge.  The  only  difference  between  the  two 
codes  is  a  charge  identification  entry  (NCHRG)  and  the 
charge  weight  and  constants  specified  for  the  conventional 
charge  instead  of  the  pressure-time  history  (WCHRG-THEXP) . 


FILEi  USL46 


DATA 


A1 


1599  16  2  9  2  2  1  /  NTIME.NSK1P«NCHRG,NQUAD,NFINE,IC0UPLE/NSUBS 
5.0E-05  -640.0  0.0  0.0  /  OELT.XLOAO.RLOAO.SURCUT 
9  /  NSHAPE 

0.0  240.0  0.0002  320.0  0.00053  350.0  0.00243  305.0  0.00294  160.0 
0.0034  115.0  0.0044  70.0  0.0054  46.0  0.00565  0.0  /  (TC(J),PC(J),Jsl, NSHAPE) 
0.0  /  DECAY 
0100000000 

0010000000/  OUTPUT  FLAGS  1-20 
4  /  NPTSHL 

1  1  2  11  2  46  3  12  /  LBOSE6.LBOSPT 
1  /  NPTSUB 


FILEi  USL47  DATA  A1 


1599  16  1  9  2  2  1  /  NTiriE,NSKIP,NCHRG,NQUAD,NFINE,KOUPLE,NSUBS 
5.0E-05  -840.0  0.0  0.0  /  OELT,XLOAD,RLOAD,SURCUT 

352.0  3.8354314E+5  1.144  3.03131E-5  -.247  /  WCHRG,PZMLT,PZEXP,THMLT,THEXP 
0100  0  0  0000 

0010000000/  OUTPUT  FLAGS  1-20 
4  /  NPTSHL 

1  1  2  11  2  46  3  12  /  LBOSEG.LBOSPT 
1  /  NPTSUB 


'i'vVVv 


7.  POSLOB  INPUT  DATA 

The  PUSLOB  code  produces  TEKTRONIX  plots  of  the 
velocity-time  histories  for  specified  points  on  the  shell  or 
substructure,  it  also  generates  velocity  punch  card  files 
which  are  converted  to  a  format  for  EASYPLOT  using  the 
PLOTCONVl  code. 


<fN) 


FILE:  PUSL46  DATA 


AX 


1599  16  I  I  1  /  MTIME,NSKIP,MSUBS,NTEK,NCARD 
5.0E-05  1000.0  1.0  /  DELT,XMULT,YMULT 
VELOCITY  PROFILE 
TIME  (MSEC) 

VELOCITY  (IN/SEC) 

NPTSHL-W 
-FWD  1/01 
11-11  0.0 
V-FWD  3/12 

3  12  -1  1  0.0  /  LBOSEG,LBOSPT,NSPHC,NANG,ANGDEG 
13  NPTSHL-U 
V-FWD  1/01 
1  1  1  1  0.0 
V-FWD  2/02 
2  2  0  1  0.0 
V-FWD  2/07 
2  7  0  1  0.0 
V-FWD  2/11 
2  11  0  1  0.0 
V-FWD  2/15 
•  2  15  0  1  0.0 

V-FWD  2/20 
2  20  0  1  0.0 
V-FWD  2/24 
2  24  0  1  0.0 
V-FWD  2/28 
2  28  0  1  0.0 
V-FWD  2/33 
2  33  0  1  0.0 
V-FWD  2/37 
2  37  0  I  0.0 
V-FWD  2/41 
2  41  0  1  0.0 
V-FWD  2/46 

2  46  0  1  0.0 
V-FWD  3/12 

3  12  1  1  0.0  /  LBOSEG,LBOSPT,NSPHC,NANG,ANGDEG 

2  NPTSHL-V 

V-FWD  l/Ol 
11-11  0.0 
V-FWD  3/12 

3  12  -1  1  0.0  /  LBOSEG,LSOSPT,NSPHC,NANG,ANGDEG 
1  NPTSUB 

PLT  CNTR 
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